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The most striking characteristic of crystalline solids is their periodicity. As a result of this feature, 
theoretical descriptions of physical phenomena in such systems are usually <iiven in wave number or 
momentum space. The reciprocal lattice of a crystal and the Fermi surface of a metal are examples. In 
a disordered system, on the other hand, there is no such periodicity and momentum space descriptions 
are much less natural. However, in such systems, physical conditions near a point r, in coordinate 
space, become independent of the conditions at a distant point r', provided that |r'-r| is Iar<je compared 
to either a characteristic mean free path or some other appropriate length. This suggests that one can 
analyze a macroscopic disordered system by averapn<i over the properties of microscopic neighbor- 
hoods. 

In the present paper we report some details of such a program. Although the point of view is of 
quite general applicability we have, for the sake of definiteness, studied so far only one type of system: 
Noninteractin^ electrons moving in the field of interacting, disordered scattering centers. We have 
focused especially on the electronic density of states. The macroscopic system is represented by an 
average over small neighborhoods. If one did not take special precautions, one would encounter one 
class of errors of the order of d/L where L is a characteristic dimension of the neighborhood, and d is a 
characteristic atomic dimension; and another class of errors of the order of 1/JV where /V is the number 
of ions. Both are too large to be tolerable for practical purposes. However, by an appropriate treatment 
of the statistical mechanics of the scatterers and by periodic repetition of the small neighborhoods, 
these errors can be avoided. The remaining errors are exponentially small in the ratio y(LIH) where y is 
of order unity and R is the smaller of the electronic mean free path or the deBroglie wavelength of the 
electrons. This exponential convergence of the small neighborhood theory promises to make it a useful 
practical method for the study of disordered systems, especially very highly disordered ones. 

Numerical examples are presented and discussed. 

Key words: Binary alloys; density of states; disordered systems; periodically continued neighbor- 
hood. 



1 . Introduction teristics at a point r depend significantly only on condi- 
The physical properties of strongly interacting disor- tions inside a rather sma11 neighborhood |r'-r|<R, 
dered systems are in general difficult to calculate since where R is of the order of a mean free P ath l or a 
simplifying symmetries, present in crystalline mated- characteristic thermal deBroglie wavelength, A,. The 
als, are not present. In this paper we present and effects of the more distant environment fall off ex- 
develop to some small extent a viewpoint which ap- ponentially with distance. This suggests that the pro- 
pears to provide a useful line of attack on the theory of perties of an infinite disordered system could be accu- 
disordered systems. rately calculated by suitable averaging over an ensem- 
We shall demonstrate and utilize the rather plausible ble of small neighborhoods. This scheme, furthermore, 
fact that in a disordered system the physical charac- should be most successful for strongly disordered 

systems in which the mean free path is short. 



In the following sections we sketch an application of 
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of states of noninteracting electrons moving in the 
static potential of disordered scattering centers. 3 

2. Demonstration of the Locality Principle 

In this section we shall sketch the demonstration that 
the physical properties near a point r of a system of 
noninteracting electrons moving in a disordered exter- 
nal potential are nearly independent of the circum- 
stances outside a characteristic range R. More specifi- 
cally, the effect of a perturbation at a distant point r' on 
physical properties at r falls off exponentially, on a 
scale given by a range of influence, R, which is of the 
order of either the mean free path or, at high tempera- 
tures, the thermal wavelength of the electrons, 
whichever is the smaller. 

We consider a very large system of volume Cl 
described by the Hamiltonian 



H = Hq+V+v; 



(2.1) 



here Ho is a periodic Hamiltonian, which for simplicity 
we take to be the kinetic energy 



tf,» = P 2 ; 



(2.2) 



V is a potential produced by disordered scattering cen- 
ters, located at the points r a , 



V=^V a {v-v a 



(2.3) 



and v is an additional small perturbation localized at a 
point r'. We shall study the effect of v on physical con- 
ditions at the point r. 

First, we shall consider the one particle Green's func- 
tion 

1 



G(r,r';£) = r 



E-H 



(2.4) 



In terms of the eigenfunctions t/r n of H, and the eigen- 
values E n , 

G(T.r';E)=2 *!(T)MT') . 

n Z-En 

Of special interest is the contracted function G(r, r; E) 
which determines the density of states by means of the 
following relation, 



n(E) 



->{ 



drG(r, r;£ + £0). (2.6) 



• ! Among fairly recent papers we mention the following: H. Schmidt, Phys. Rev. 105,425 
(1957): S. F. Edwards, Phil. Mag. 3, 1020 (1958); ibid. 6, 617 (1961); Proc. Roy. Soc. A267, 
518 (1962): J. L, Beeby, Proc. Roy. Soc, A279, 82 (1964); Phys. Rev. 135A, 130 (1966): P. 
Seven. Phys. Rev. 151, 539 (1966); ibid., 156, 809 (1967). 



We see that we may regard the unintegrated quantity 






n(r,E) s-i ImC(r, r; £ + $) 
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■ 

(2.7) - 



as the local density of states density. 

We would like to estimate the effect of introducing 
the weak additional potential v, at the point r', on n(r, 
E) at the point r. The quantity of physical interest is 
then 



— r -rrGfr, r;E + i0) 
bv(r ) i 

= <G(r, r';£+i0)G(r\r;£+*0)%=o (2.8) t] 

where the brackets () denote configuration average. 
The equation (2.9) follows directly from the equation of 
motion of G. 

Now it is well known that, for weak random poten- 
tials V a , 

i 






<G(r,r';E+t0)> 



477|r— r'l 



,iV£|l— r'U-|r- r'|/21(«), 



(2.9) 



where 1(E) is the mean free path of an electron of ener- 
gy E. Similarly, one can show from (2.8) that for large 



Mr') 



G(r,r;E+i0))\-~e 



r'|//(£) 



(2.10) 



Thus we see that when |r'-r| > 1(E), a change of poten- 
tial at r' has a negligible effect on the quantity n(r, E). 
Next we show a similar locality effect, this time not 
due to a finite mean free path but due to elevated tem- 
perature. We consider a system of independent elec- 
trons and use Boltzmann statistics for simplicity. We 
take as Hamiltonian 



H = H + v, 



(2.11) 



where, in comparison with (2.1), we have eliminated the 
disordered scattering potentials. We write the partition 
function as 



Z=JdrZ(r), 
Z(r) = (r|e-**|r). 



(2.12) 



(2.13) 



Our interest now is in the influence of a small perturba- 
tion v at r' on Z(r). By standard perturbation theory one 
can show that, for large | r ' — r | , 



77' 



5/2 



8Z(r) 

dv(r') 16£ 3 / 2 |r- 



-|r-r'| 2 /A7 



(2.14) 



444 



where \i is the thermal wavelength 

\ t = pH*(=hl(2mkT)U*). 



(2.15) 



Of course at high temperature and with disorder, the 
characteristic "range of influence" R will he either a 
representative mean free path Z, or A*, whichever is the 
smaller. 

3. The Periodically Continued Neighborhood 

We shall now develop a concrete method of calcula- 
tion, based on the locality principle of the previous sec- 
tion. We shall concentrate on the density of states n(E). 

The most straightforward way would be as follows. 
We imagine the large disordered system as given. We 
choose at random a large number of points R/ and sur- 
round them by spheres of radius p considerably larger 
than the "range of influence" R, but not too large. Each 
sphere is surrounded by an infinite wall (see fig. 1). The 
electronic Hamiltonian for each sphere is given by 



where 



H=H Q +V 



V=\XV a (r-r„) |r-R,|^p 



+ < 



r>p 



(3.1) 



(3.2) 



We then calculate for each sphere the density of 
states density at its center, n(Ri,E). Because of the rela- 
tively small size of the neighborhood this is a much 
more manageable problem than n(E) for the macro- 
scopic disordered system. In view of the locality princi- 
ple, n(¥LiJZ) is only insignificantly affected by the 
presence of the infinite wall. Hence, the density of 




states for the macroscopic system is given by 
n(E) - (ln(Ri£) 9 



(3.3) 



where the bar denotes an average over many R/. 4 

The practical drawbacks, for finite p/Ro, are two: (1) 
The replacement, in the outside region, of the actual 
potential, 2 V a , by an infinite repulsive barrier is a quite 
drastic change and unless p/R is very large, will cause 
sizeable errors. (2) This method does of course not lead 
to the exact results in the special case of a vanishing or 
periodic potential. 

These drawbacks can be largely overcome by using, 
instead of a finite and bounded neighborhood, as in 
figure 1, a periodically continued finite neighborhood 
as indicated in figure 2. We choose a fundamental cell, 
say the cube, of volume 11/, = L 3 , 

L 



WU 



L>2R 



(3.4) 



and construct the space-lattice generated by it. Let us 
call t {v) the lattice translation vectors. Then we place 
any number of scatterers in some definite configuration 
c in the fundamental cell and populate the other cells in 
the identical way (see fig. 2). The electrons now move 
in a periodic cubic lattice. Their energy eigenvalues 
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Figure 1. Spherical neighborhood centered at R/. 



Figure 2. Periodically continued neighborhood. 



4 If we would use the total densities of states of the spherical neighborhoods, rather than 
the quantities n(Ri,E) computed at their centers, we would incur large emus behaving as 
L~\ due to the presence of the infinite wall boundaries. 
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are, because of the relatively small size of the unit 
cell, much more amenable to calculation than the 
eigenvalues of the macroscopic disordered system. Let 
us call the density of states, corresponding to the con- 
figuration c and total volume ft, n c (E). Then our approx- 
imation for the density of states of the actual macro- 
scopic system will be 



n(E)~ ^w c n c (E\ 



(3.5) 



where w c are weights, which we shall presently discuss, 
and the approximate equality, ~, will signify accuracy 
to within terms exponentially small in L/2R. 5 

We shall now describe a suitable choice of w c and 
later demonstrate that it leads to the claimed accuracy. 
Let us suppose that in the actual, macroscopic system 
under consideration there are given two-body forces 
between the scatterers, 



<Pap = <Pap(r a — rp), 



(3.6) 



which vanish beyond a range a which is much smaller 
than L. Let us suppose further than the scatterers obey 
classical statistics. Thus, in the macroscopic system 
the probability of a given configuration c = fri, r 2 , . . .) 
is given by the grand canonical weight function, 



v c ^ = A icc) exp J —/3 



« 2 <P«0-]£ fJ-jNj 



(3.7) 



a^p 



where A {oo) is the normalization constant, /Xj is the 
chemical potential of scatterer of type j, and Nj is the 
total number of scatterers of this type. 6 

A suitable choice of w c , for the periodically continued 
neighborhood, is obtained as follows. We take the unit 
cell and associate with each original position vector r, 
the infinite set of vectors 



p(*0 = r + t (p) 



(3.8) 



This corresponds to the topology of a three dimensional 
torus and is illustrated in figure 3 in one dimension. 
The weight w c is then determined by the equation 



w c 



A exp \ — /3 



i|)v^(rctf)-2>ty]) (3.9) 



5 Provided the forces between the scatterers are sufficiently short range. 

6 For the macroscopic system, one could of course equally well use a canonical distribu- 
tion. However, in our local neighborhood theory, this would lead to unacceptable errors of 
order 1//V, where Visa mean number of scatterers in fl/.. 




Figure 3. Schematic representation of the toroidal topology. The 
circumference is L. 

where r a p is the shortest vector contained in the set 
r a (/x) — r/3 { ^\ For example, r 12 is shown schematically 
in figure 3. The normalization constant A is chosen so 
that 



w c 



i, 



(3.10) 



where the sum goes over all configurations, including 
all possible numbers of the scatterers. 

It can be shown that, with the choice (3.10), the cor- 
relation functions, ft s (ri .... r s |L), in the finite toroidal 
system, differ negligibly from those of the infinite 
system, n s (r u . . . . r s | °°), for values of |r a — r^l^L/2 
and s up to s^L/a. The error of a given correlation func- 
tion behaves as exp [— a s L/a], where a s is of order unity. 
This fact assures that, if L > a, the statistical distribu- 
tions of the ions in any neighborhood of size < L/2 are 
practically identical in the ensemble of periodically 
continued neighborhoods and in the macroscopic 
system. Hence, in view of the locality principle demon- 
strated in section 2, the density of states n(E) of the in- 
finite system may be determined, via eq (3.5), from the 
density of states in the periodically continued neighbor- 
hoods. The error will be exponentially small in the 
quantity L/R or L/a, whichever is the smaller. Typically 
both a and R are of the order of 1 A, so that one must 
work with neighborhoods of dimensions of several in 
order to obtain quantitatively useful results. 

It is evident that this method will give exact results 
for perfectly periodic systems. Since it is also very ac- 
curate for highly disordered systems (small R), it should 
give good answers for most intermediate situations. 

The ensemble of systems defined by eq (3.9) is a 
grand canonical ensemble in which the volume £Il(=L 3 ) 
is fixed while the numbers of particles Nj assumes all 
possible values. For a single species of atoms it has 
been found practically preferable to work with a 
periodically continued isothermal-isobaric ensemble, 
in which TV is fixed but the volume ft/, is variable. Here 



446 



a unit cell of volume Hl, with the TV atoms in a configu- 
ration c must be given the weight 

^ nL * c= nI exp l ^[1 2 *+(*+) -R&l }(3.ii) 

where P is the pressure. The density of states of the ac- 
tual macroscopic system is approximated by 

n(E) - [ d£li. 2 wa L ,ME; 11/,; c) (3.12) 

J c 

where n(E; CIl; c) is the density of states, per unit cell 
flz,, in the periodically continued neighborhoods. The 
normalization A of wsi L , c 

[ dCl L %wn L , c =N tot IN (3.13) 

j r 

where Ntot is the total number of scatterers in the 
macroscopic system. Again the convergence of (3.13) to 
the exact result is exponential, as for the case of the 
grand canonical ensemble. 

4. Numerical Illustration 

To illustrate the theory of the previous sections we 
have numerically studied the following model of a one- 
dimensional alloy: 



The density of states curves was calculated for each 
of the 2 6 configurations and then averaged. For L — 10 
a sampling procedure was used. The results are shown 
in figure 5, together with the exact result obtained for 
the infinite system using the Schmidt method. 3 Our ap- 
proximation reproduces quite accurately all the details 
of the density of states structure, even for L = 6. To ob- 
tain similar accuracy from a single randomly populated 
chain would require a length of the order of L ~ 10 3 . 



H = --j-r l + ^KMx-ad) 



(4.1) 



where the distance between potentials, d, was taken as 
1, and k a was taken, with equal probability, as —2 and 
—4. (This corresponds to a high temperature limit for the 
two different kinds of scatterers, an idealized model for 
- crystalline Cu Au well above the ordering temperature.) 
Here the grand canonical method was appropriate. Two 
calculations were performed, with L = 6 and L=10. Let 
us take the case L = 6. A typical, periodically continued 
neighborhood is shown schematically in figure 4. Each 
configuration of potentials in the fundamental interval 
< x ^ 6 has the same weight, 2~ 6 . 

In this simple example all correlation functions are 
evidently exact, as long as the crystal sites considered 
are all contained in a single unit cell. 
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FIGURE 4. A typical periodically continued neighborhood of the 
illustrative model. 



FIGURE 5. Integrated density of states for binary alloy. The solid line 
is the exact result. 



5. Concluding Remarks 

It is commonplace to emphasize the theoretical dif- 
ficulties caused by disorder. In the present paper we 
draw attention to a favorable feature: For a highly disor- 
dered system the physical properties near a given point 
depend only on circumstances in a small neighborhood, 
whose dimension is of the order a mean free path. Con- 
sequently, a macro-system may be treated, so to speak, 
neighborhood by neighborhood and the macro-problem 
can be reduced to an ensemble of micro-problems. We 
have shown that the errors of such a procedure can be 
made to vanish exponentially with the size of the 
neighborhood. A numerical model calculation bears out 
these considerations. 

A more complete account of this program will be 
published elsewhere. 

(Paper 74A3 -610) 
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